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Abstract 

We apply a numerical minimum action method derived from the Wentzell-Freidlin theory of large 
deviations to the Kardar-Parisi-Zhang equation for a growing interface. In one dimension we find 
that the switching scenario is determined by the nucleation and subsequent propagation of facets 
or steps, corresponding to moving domain walls or growth modes in the underlying noise driven 
Burgers equation. The transition scenario is in accordance with recent analytical studies of the 
one dimensional Kardar-Parisi-Zhang equation in the asymptotic weak noise limit. We also briefly 
discuss transitions in two dimensions. 
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I. INTRODUCTION 



The large majority of natural phenomena are characterized by being out of equilibrium. 
This class includes turbulence in fluids, interface and growth problems, chemical reactions, 
processes in glasses and amorphous systems, biological processes, and even aspects of eco- 
nomical and sociological structures l2i- In this context there is a continuing interest in 
the strong coupling aspects of stochastically driven non equilibrium model systems 0, 4|. 
Here the dynamics of complex systems driven by weak noise, corresponding to rare events, 
is of particular interest. The issue of different time scales characterizes many interesting 
processes in nature. For instance, in the case of chemical reactions the reaction time is often 

n 

orders of magnitude larger than the molecular vibration periods [5| . The time scale separa- 
tion problem is also encountered in the case of conformational changes of biomolecules , 
nucleation events during phase transitions, switching of the magnetization in magnetic ma- 
terials 71, 81 , and even in the case of comets exhibiting rapid transitions between heliocentric 
orbits around Jupiter j^. 

In the weak noise limit the standard Monte Carlo method or direct simulation of the 
Langevin equation becomes impractical owing to the large separation of time scales and 
alternative methods have been developed. The most notable analytical approach is the 
formulation due to Freidlin and Wentzel which yields the transition probabilities in terms of 



an action functional [lO|]. This _approach is the analogue of the variationa' 



by Machlup and Onsager [ll, ll2|, see also work by Graham et al. 



13 



principle proposed 
14l ] and Dvkma n |l5 |. 



The Freidlin- Wentzel approach is also equivalent to the Martin-Siggia-Rose method 16|] in 



the weak noise limit of the path integral formulation [17|, |18|, [191, |20|, |2l| . In order to overcome 
the time scale gap various numerical methods have also been proposed. We mention here 

3, 25, 26]. 
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the transition path sampling method [22| and the string method 

A particularly interesting non equilibrium problem of relevance in the nanophysics of 
magnetic switches is the influence of thermal noise on two-level systems with spatial degrees 



of 



reedom 



m 



271]. In a recent paper by E, Ren, and Vanden-Eijnden |28|], see also Ref. 
291 ]. this problem has been addressed using the Ginzburg-Landau (GL) equation driven by 



thermal noise. Applying the field theoretic version of the Onsager-Machlup functional 



11 



121 ] in the Freidlin- Wentzell formulation [10|] , these authors developed the so called minimum 



action method in which they implemented a powerful numerical optimization techniques for 



2 



the determination of the space-time configuration which minimizes the Freidhn-Wentzell 
action. The minimizers correspond to the kinetic pathways and the associated action yields 
the switching probabihties in the long time-low temperature limit. In the picture emerging 
from the numerical study the switching between metastable states is due to noise induced 
nucleation and subsequent propagation of domain walls across the sample. Subsequently, we 
supplemented the work by E et al. and presented a dynamical description and analysis of 
the non equilibrium transitions in the noisy ID GL equation for an extensive system based 
on a weak noise canonical phase space formulation of the Freidlin-Wentzel or Martin-Siggia- 



Rose methods 
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31| . We derived propagating nonlinear domain wall or soliton solutions 



28 



29|. 



of the resulting canonical field equations with superimposed diffusive modes. The transition 
pathways are characterized by the nucleations and subsequent propagation of domain walls. 
We discussed the general switching scenario in terms of a dilute gas of propagating domain 
walls and evaluated the Arrhenius factor in terms of the associated action. In conclusion we 
found excellent agreement with the numerical studies by E et al. 

The noise-driven GL equation belongs to the class of so-called gradient systems where the 
drift term in the Langevin equation can be derived from a free energy functional. Regarding 
the kinetic transitions this property implies the existence of an underlying free energy land- 
scape in which the optimal pathway proceeds via saddle points, yielding the corresponding 
Arrhenius factor. We note that this interpretation implies a fluctuation-dissipation theorem 
relating the strength of the noise to the kinetic transport coefficient. There is, however, 
another interesting class of stochastic model systems characterized by Langevin equations, 
where the drift term cannot be associated with a free energy functional. Those are the so 
called non-gradient systems for which the interpretation of pathways in a free energy land- 
scape fails and has to be replaced by pathways in an "action landscape". In recent work. 



see e.g. Refs. 



32 



33], where earlier references can be found, we have addressed a non equi- 



librium model falling in the class of non-gradient systems, namely the Kardar-Parisi- Zhang 
(KPZ) equation or, equivalently, noisy Burgers equation describing, for example, a growing 
interface. Using the weak noise canonical phase space method alluded to above, we find that 
the kinetic pathways correspond to nucleation and propagation of localized growth modes 
with superimposed diffusive modes. The growth modes together with the diffusive modes 
carry an action, yielding the transition probabilities. 

The purpose of the present paper is to attempt to substantiate the weak noise growth 
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mode approach to the KPZ equation by a direct numerical optimization employing the 



minimum action method developed by E, Ren, and Vanden-Eijnden 
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281], see also Refs. 



351 ] ■ Similar to the GL case we find in ID that the switching scenario is determined by 



the nucleation and propagation of growth modes. We are also able to account numerically for 
the associated transition probabilities. The paper is organized in the following manner. In 
Sec. [TTlwe briefiy review the KPZ equation and the analytical results obtained by the weak 
noise canonical phase space approach. In Sec. Illll we introduce the minimum action method 
and establish the connection with the phase space method and path integral formulations. 
In Sec. llVl we discuss the numerical implementation of the Freidlin-Wentzel scheme. In Sec. 
IVl we present the numerical results for various switching scenarios. In Sec. |Vl] we offer a 
heuristic discussion of the numerical data based on the analytical phase space method. In 
Sec. IVIII we briefly discuss kinetic pathways in 2D and present some switching scenarios. 
Sec. IVIIIl is devoted to a summary and a conclusion. 



II. THE KPZ EQUATION 

In this section we review the KPZ equation and apply the weak noise method. The KPZ 
equation describes an intrinsic non equilibrium problem and plays in some sense the same 
role as the Ginzburg-Landau functional in equilibrium physics |2|, |36|. The KPZ equation 



was introduced in 1986 in a seminal paper by Kardar, Parisi, and Zhang 



see also Ref. 



38l . l39| , and purports to describe non equilibrium aspects of a growing interface [4, |40| . In 
a Monge representation the KPZ equation for the stochastic time evolution of the height 
fleld h{r, t) has the form 

^ = uV'h + ^{Vhr-F + v, (2.1) 
{Vv){r,t) = A6''{r)6{t). (2.2) 

Here the damping coefficient or viscosity u characterizes the linear diffusion term uV'^h, the 
growth parameter A controls the strength of the nonlinear growth term (A/2)(V/;,)^, F is 
the imposed drift term and is a constant here, and r] a locally correlated white Gaussian 
noise, modelling the stochastic nature of the drive or environment; the noise correlations are 
characterized by the noise strength A. 



4 



A. Burgers and Cole-Hopf equations 



In the growth mode analysis of the KPZ equation the local slope of the growing interface 
given by the vector field, 



u = Vh, 



(2.3) 



is of importance. In terms of u the 
driven by conserved noise 4l|, |42 



43 



iPZ equation then maps onto the Burgers equation 



44| 



du 

'dt 



z/V^u + A(u ■ V)u + Vf]. 



(2.4) 



In the deterministic case for r ? = the Burgers equation has been used to study irrotational 



fluid motion and turbulence 



45 



46 



47 



49] ; it has also played a role in astrophysics in 



the context of large scale structures in the universe 50, l5l[ |. 



Another quantity of importance in our analysis of the KPZ equation is the diffusive field 



w defined by the non-linear Cole-Hopf transformation 38 

exp 



52 



53| 



w 



A, 

2u 



(2.5) 



In terms of w the KPZ equation maps onto a linear diffusion equation driven by a multi- 
plicative noise, here denoted the Cole-Hopf (CH) equation, 

dw 
'dt 



^2 A ^ 

uV w wr H wri. 

2z/ 2z/ 



(2.6) 



In the absence of noise for rj = the CH equation reduces to the linear diffusion equation 
and is readily analyzed permit ting a complete discussion of the KPZ and Burgers equations 



in the deterministic case 
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38| . In the noisy case a path integral representation maps 



the CH equation and consequently the KPZ equation onto a model of a directed polymer 
in a quenched random potential. The disordered directed polymer constitutes a toy model 
within the spin glass literature and has been analyzed by means of replica, Bethe ansatz. 



and functional renormalization group techniques 
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54 



55 



56|. 



B. Scaling properties 

Most work on the KPZ equation has addressed the scaling issues. For completeness we 
summarize the salient features here. The KPZ equation lives at a critical point and conforms 



to the dynamical scaling hypothesis [4, [571, l58|, l59|] which in terms of the height correlations 
assumes the form: 



{hh)ir,t) = r'^fit/r^). (2.7) 

Here ( is the roughness exponent, z the dynamical exponent, and / the associated scaling 
function. The exponent C is a measure of the roughness of the interface, e.g., for C = the 
interface is flat, for ( = 1/2 the interface exhibits a random walk profile, {hh){r) ~ r. The 
exponent z is a measure of the dynamical scaling, e.g., for diffusive behavior z locks onto 2; 
for a ID growing interface z = 3/2. 

In order to extract scaling properties the initial analysis of the KPZ equation was based 



on the dynamic renormalization group (DRG) method, previousl y d eve. 



oped and applied to 



dynamical critical phenomena and noise driven hydrodynamics j41|, |42, |60|| . An expansion 
in powers of A in combination with a momentum shell integration yields to leading order in 
d — 2 the DRG equation dg/dl = P{g), with beta-function P{g) = {2 — d)g + const. g"^ [Gil. |62|. 



Here g = A^A/z/^ is the effective coupling strength and / the logarithmic scale parameter. 
Below the lower critical dimension d = 2 the DRG flow is towards a strong coupling fixed 
point with scaling exponents ( = 1/2 and z = 3/2 in d = 1. Above d = 2 a kinetic phase 
transition line delimits a strong coupling regime from a weak coupling regime. In the strong 
coupling regime the DRG fiow is towards a still poorly understood strong coupling fixed 
point with unknown scaling exponents and scaling function. In the weak coupling regime 
the DRG flow is towards a weak coupling flxed point with scaling exponents z = 2 and 
C = {2 — d)/2. On the transition line z = 2 and C = 0? see e.g. Ref. [33|. In Fig. [T]we have 
depicted the scaling properties in a plot of the coupling strength g versus the dimension d. 

We note two further properties of the KPZ equation which are also relevant in a scaling 
context. Firstly, subject to a Galilean transformation the equation is invariant provided we 
add a constant slope to the height field h and adjust the drift F accordingly, i.e., 

r ^ r - Au°t, (2.8) 
h^h + u° -r, (2.9) 
F^F+(A/2)u°-u°; (2.10) 

note that the slope field u and diffusive field w transform according to u ^ u + u° and 
w w exp[[X/2h')u^ ■ r], respectively. From a simple scaling argument and also following 



from the DRG analysis the Galilean invariance implies the scaling law 



(2.11) 



relating the roughness and dynamic scaling exponents. The Galilean invariance is a funda- 
mental dynamical symmetry specific to the KPZ equation, delimiting the KPZ universality 
class. Secondly, a fluctuation dissipation theorem is operational in ID in the sense that the 



stationary 



solution 
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bkker-Planck equation associated with the KPZ equation admits the explicit 
6J 



Po{h) oc exp 



'a 



dx{VKY 



[212) 



The Gaussian form of the distribution shows that the slope u = Vh fluctuations are uncor- 
related and that the height field h = udx' performs a random walk in x. Note also that 
the distribution is independent of the non- linear growth parameter A. 



C. Weak noise method 



Whereas the DRG approach, based on an asymptotic expansion about the critical di- 
mension d = 2, deals with the long time - large distance scaling properties of the KPZ 
equation, the asymptotic weak noise approach addresses the stochastic growth morphology 
or many body aspects. The weak noise or canonical phase space method focuses on the 
noise strength A as the relevant parameter in the problem. In the absence of noise for rj = 
or A = the morphology of the deterministic KPZ equation decays subject to a transient 
pattern formation; in ID corresponding to cusps connected by parabolic segments [38|. In 
the presence of even weak noise the KPZ equation is eventually driven into a stationary 
stochastic state; the cross-over time diverging in the limit of vanishing noise. In this sense 
the noise strength A is a singular parameter and a weak noise approach asymptotic in A. 

The weak noise canonical phase space approach is implemented by applying an eikonal or 
WKB approximation to the Fokker-Planck equation associated with the Langevin equation. 
Viewing the Fokker-Planck equation as an imaginary time Schodinger equation the scheme is 
equivalent to the well-known WKB or semi-classical approximation in quantum mechanics, 
where the wave function \E' is related to the classical action 5 by \[' oc exp[iS/h], h being the 
Planck constant. In quantum mechanics the quantum fluctuations characterized by orbitals 
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are then in the correspondence hmit ^ ^ replaced by orbits as solutions to the classical 
equations of motion following from the action 5*. 

In the weak noise approach the point of departure is a general Langevin equation of the 
form 



^ = -F(x)+r7(t), 



(2.13) 
(2.14) 



for sim plic ity we consider a single random variable x{t); for the more general case see e.g. 



Refs. 33|, l6J]. Here F{x) is a general non-linear drift, and rj an additive white noise 
correlated with strength A. In order to implement the weak noise approximation we consider 
the equivalent Fokker-Planck equation for the distribution P{x,t): 

. dP 



dt 



(2.15) 



Interpreting Ad/dx as a momentum operator, P as an effective wave function, and A as 
an effective Planck constant, Eq. (12.151) has the form of an imaginary time Schrodinger 
equation. Consequently, in the weak noise limit it is suggestive to introduce the WKB or 
eikonal approximation 65| 

s{x,Ty 



P{x, T) (X exp 



A 



(2.16) 



To leading order in A the action S then obeys a principle of least action ^S* = as expressed 
by the Hamilton- Jacobi equation dS/dt + H{x,p) = with associated canonical momentum 
p = dS/dx 66|, l67|]. The Hamiltonian (energy) takes the form 



H 



]^p^ -pF = ^p{p-2F), 



(2.17) 



yielding the coupled Hamilton equations of motion 



dx 

'dt 

dp 



-F + p, 
dF 



dt ^ dx 



(2.18) 
(2.19) 



Finally, the action associated with an orbit from Xi to x in the transition time T is given by 



S{x,T) 



px,T 


dx 


I dt 


^Tt - ^. 





(2.20) 



or inserting the equations of motion for x 

S{x,T) = \ dtp\ (2.21) 

The issue of solving the stochastic Langevin equation (12.1 3p or, equivalently, the determin- 
istic Fokker-Planck equation (I2.15P in the weak noise hmit A ^ is then replaced by, as 
first step, solving the coupled equations of motion (I2.18H2.191) for an orbit from an initial 
configuration Xj at time t = to a final configuration x at time t = T. In the next step we 
evaluate the action S associated with the orbit and infer from the WKB ansatz (I2.16P the 
transition probability for the specific transition. We note that the noise in Eq. (12.131) has 
been replaced by the canonical momentum p and that p is a dependent variable which has 
to be chosen in accordance with the initial and final values of x and the imposed transition 
time T. 

In a phase space representation the zero-energy manifolds p = and p = 2F intersecting 
at a hyperbolic saddle point play an important role in determining the long time stationary 
distribution Po{x) = limT^oo P{x,T). Initially an orbit from Xi to x moves along the the 
transient zero-energy manifold p = towards the saddle point. This part of the orbit repre- 
sents the transient motion. As time progresses the orbits bends away from the saddle point 
and is attracted to the stationary submanifold p = 2F. This part of the orbit corresponds to 
the cross-over to a stationary random motion. In the limit of a long transition time the orbit 
from Xi to X passes close to the saddle point, at which the large waiting time ensures the 
Markoff property. In Fig. [2] we have sketched the {x,p} phase space showing the zero-energy 
submanifolds, the saddle point and an orbit from Xi to x in transition time T. 



D. Growth modes 



In the KPZ case the weak noise scheme is most easily implemented for the CH equation 
(12. 6p driven by multi plic ative noise. This requires an extension of the weak noise method 



discussed in Ref. 



32 



33| . Introducing the wavenumber parameters 



k = {\F/2v^fl\ 



(2.22) 
(2.23) 
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setting inverse length scales, we find the weak noise Hamiltonian 

H= j (fx {p[vV'^ -vk:^]w+{l/2)kl{wpf) , (2.24) 



and associated field equations 

dxv 

— = u[V^w - k'^w] + klw'^p, (2.25) 

^ = -r,[v'p-ep]-kyw, (2.26) 

determining orbits in a {w,p} phase space. Likewise, one infers the action 



1, 



S{w,T) = ^klJ d'^xdt{wpf, (2.27) 

yielding the transition probability to leading asymptotic order in A 

P(^,T) ocexpL^^^l; (2.28) 

note that on the p = manifold Eq. fl2.25p reduces to the deterministic CH equation for 
7] = 0. 

The equations of motion fl2.25H2.26l) serve two purposes. On the one hand, a solution or 
orbit in phase space from an initial configuration Wi{v) at time t = to a final configuration 
w{y) at time t = T with p as an adjusted noise field yields an action S and thus a contribution 
to the transition probability P{w, T). Secondly, the solution w(r, t) interpreted as a classical 
orbit also provides a growth morphology for the CH equation. The deterministic growth 
or evolution of the diffusive field w then corresponds to a growth morphology for the KPZ 
equation by means of the inverse Cole-Hopf transformation 

h={l/ko)\nw. (2.29) 

Likewise, the transition probability P{h,T) is given by 

P(h, T)= fjl dwS{h - (l/ko) logw)P{w, T). (2.30) 

r 

The growth morphology follows from the coupled nonlinear field equations (12.25112.2^ . Ow- 



ing to the negative diffusion coefficient the equations are numerically unstable, see Ref . [68| , 
however, searching for localized instanton or soliton type solutions we note that on the p = 
and p = uw submanifolds the static equations reduce to the static diffusion equation and 
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the static nonlinear Schrodinger equation, well known in the context of dark solitons in Bose 



condensed atomic gasses 



68| 



V w = k 



(2.31) 
(2.32) 



E. Domain walls in ID 



In ID Eqs. (12.31112.3^ admit the static solutions w± oc cosh^^ kx for the diffusive field 
w. These modes correspond to cusps in the height field, h± = ±{l/ko) ln(cosh A;x), and to 
static domain walls or solitons in the local slope field: 



k 

u±{x) = ±— tanh/cx. 
ko 



(2.33) 



The right hand domain wall, u^{x) = (k/ko) tsmhkx, is associated with the p = manifold 
and carries zero energy and zero action. This mode is the well-known viscosity-smoothed 
shock wave solution of the static noiseless Burgers equation uV'^u + XuVu = 0, as easily 
seen by inspection 4J]. The left hand domain wall, U-{x) = — (/e/Zeq) tanh lives on the 
p = uw manifold and carries a finite action 



(2.34) 



the static domain walls are depicted in Fig. [31 By applying the Galilean transformation 
(12.8112. 10() the static domain walls can be boosted to a finite propagation velocity and we 
obtain the moving domain walls or growth modes 



k 

u±{x, t) = ±— tanh k(x — Xu^t) — uq. 
ko 



(2.35) 



The propagating domain walls form the basic building blocks in the construction of a growth 
morphology. Considering a dilute gas of non overlapping growth modes of different ampli- 
tudes or "charges" ki, where a positive charge corresponds to a right hand domain wall and 
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a negative charge to a left hand domain wall, we obtain the global solution 33 1 

u{x, t) = jr 5^ tanh \ki\{x — Xi(t)), (2.36) 



1 V ^ k- 

h{x, t) = TT / , Tj~\ ln(cosh \ki\{x — Xi{t))), (2.37) 





t 





Xi{t) = / Vi{t')dt' + xl (2.38) 



Mt) = -^Ykitanh\ki\{x,it) - xiit)); (2.39) 
ko j-^ 

note that the neutrality condition fcj = ensures that the interface is flat at infinity. 
This condition, however, allows for an offset in h corresponding to propagating facets. 

The interpretation of the time dependent growth morphology is straightforward. For a 
dilute gas of growth modes the velocities adjust to constant values after a transient period 
and the growth modes move ballistically. Moreover, superimposed on the growth modes is 
a gas of diffusive modes following from a linear analysis of the equation of motion about the 
domain wall solutions, see Ref . |69j] . In Fig. H] we have depicted a three domain wall growth 
configuration composed of interconnected propagating domain walls, two right hand domain 
walls and one left hand domain wall. We also show the resulting morphology in the height 
field corresponding to moving steps or facets. 

In order to make contact with the stochastic interpretation we prepare the interface in 
a specific initial state h{x, 0) characterized by a gas of growth modes plus diffusive back- 
ground. By also assigning an appropriate noise field p(x, 0) corresponding to the nucleation 
of growth modes this configuration propagates ballistically forward in time to a specific finite 
configuration h{x,T). Only the left hand domain walls corresponding to negative charges 
carry an action. For a dilute domain wall gas, ignoring the diffusive contribution, this action 
is additive, i.e., 

ki<0 

yielding the transition probability 

S{h,T) 



P{h, T) cx exp 



(2.41) 



A 

For illustration consider the two-domain wall configuration depicted in Fig. [5l This pair 
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mode has the form 

k 

u{x, t) = —[u+{x — vt — xi) + U-{x — vt — X2)], (2.42) 

and moves according to the domain wall matching condition following from the Galilean 
symmetry with the velocity v = —Xk/k^. Since the pair mode in the slope u corresponds to 
a moving step in h the propagation across the system either subject to periodic or bouncing 
boundary condition corresponds to adding a layer to the interface; the mode thus corresponds 
to a specific growth situation. The mode moves ballistically with an action given by Eq. 
fl2.34p carried by the left hand domain wall; note that the right hand domain wall partner 
carries zero action. In time T the mode moves the distance L = vT and we obtain the 
transition probability 



P{L, T) oc exp 



4z/ L 



(2.43) 



3A2AT2 

We conclude that the step in h performs a random walk with mean square displacement 

< >oc (A2A/z/)2/3t2/^ (2.44) 
characterized by the dynamical exponent z = 3/2. This result is in accordance with estab- 



lished scaling results for the KPZ equation in ID, see e.g. Ref. 



391] . The facet in the height 



field corresponding to the pair growth mode thus performs superdiffusion [70|. 

III. THE MINIMUM ACTION METHOD 

In this section we discuss the basis for the minimum action method characterized by the 
Freidlin-Wentzel action and the connection to equivalent formulations in non equilibrium 
physics. 

A. Preidlin Wentzel scheme 

The point of departure for the Freidlin-Wentzel (FW) scheme is a generic Langevin equa- 
tion for a set of stochastic variables, {x„}, driven by additive white Gaussian noise 

^ = -F„({x™})+r^„(t), (3.1) 
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where the noise is distributed according to 



Pi{Vn)},T) oc exp 



2A 



(3.2) 



A heuristic derivation of the Freidhn-Wentzel action, the basis for the minimum action 
method, follows in the weak noise limit by simply replacing the noise ?7„ in Eq. fl3.2p by 
dxn/dt + Fn yielding 



P({x„},T) oc exp 



1 



2Aj, 

Expressing P({x„},T) in the WKB form 

P{{xn},T) OC exp 
we readily identify the Freidlin-Wentzel action 



T 



doc 



(3.3) 



^1 



FW 



S 



FW 



A 



doc f) , , 

dt + " 



(3.4) 



(3.5) 



for rigorous details see Refs. [lO|, l28 |. 

The minimum action method then corresponds to minimizing the action S'fw subject to 
an initial condition {x„(0)}, a final condition {x„(T)}, and a given transition time T. The 
method thus identifies the minimum action path in the "action landscape". The method 
works both for gradient systems where F„ can be derived from a free energy. 



(3.6) 



including e.g. the GL case and non-gradient systems like the KPZ equation. 



B. Martin Siggia Rose scheme 

The Martin-Siggia-Rose (MSR) scheme Q, 0, Q, 

ing point the Langevin equation (13. ip . For simplicity we consider, however, only a single 
stochastic variable x{t). For the transition probability P{x,T) we have by definition 



19 



20 



2l| also takes as its start- 



P{x,T) = {6{x-x{T)))„ 



(3.7) 



where we average over the noise t] driving the Langevin equation. Incorporating the Langevin 
equation determining the evolution of x{t) as a delta function constraint, averaging over the 
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noise rj according to Eq. (13.21) . noting that the change of variable from dx/dt to x yields the 
Jacobian J = exp[(l/2) J dtdF/dx], and finally setting p ^ p/A we obtain the functional 
phase space integral \7l\\ 

■ SuSR 



P{x, T) oc / Y\ dxdp 6{x — x{T)) exp 
t 

where the MSR action is given by 

>S'msr = I dt 

with MSR Hamiltonian 

-f^MSR = 



A 



(3.8) 



dx 



(3.9) 



i 2 ^ iAdF 

— p — pF -\ — . 

2^ ^ 2 dx 



(3.10) 



Since p appears quadratically it can be eliminated by a Gaussian integration and we arrive 
at the configuration space path integral 



P{x,T) (X J Y\_dxS{x — x{T)) exp 



S' 
A 



(3.11) 



with action 



S = - I dt 



dx ^ ^\ ^ dF 
dt J dx 



(3.12) 



We note that this form holds for arbitrary noise strength. The path integral is a formal 
solution of the Fokker-Planck equation. In the aymptotic weak noise limit A — only the 
saddle point in the path integral contributes. Ignoring the Jacobian contribution AdF/dx 
we recover the FW result in Eq. (13.31) in the case of one variable. 



C. Quantum analogue and phase space method 



Contact with the Fokker-Planck equation fl2.15p is easily achieved by 



noting 



thatEq. flXS]) 



73|. Introducing 



has the form of a Feynmann path integral with Planck constant A [71|, 
the momentum operator p = —iAd/dx the quantum Hamiltonian operator takes the form 

+ 

order 



H 



2 rfa;2 



■.dF\ 

lA— + 

dx y 2 dx 



(3.13) 
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where the ordering in the term {iAdF/ dx)oidev remains to be fixed. Choosing the symmetrical 
Weyl ordering {dF/dx)ordeT = {l/2){Fd/dx + dF/dx) the Schrodinger equation associated 
with H, 

dP 

tA— = HP, (3.14) 

then reduces to the Fokker-Planck equation (12.151) . Finally, formally rotating p, p ^ ip, we 
obtain a real path integral representation for P, 

S" 



P{x,T) (X / JJ^ cia;(ip5(a; — x(T)) exp 



A 



with action 

S= I dt 

and Hamiltonian 



dx 



(3.15) 



(3.16) 



In the weak noise limit A — the Jacobian contribution in Eq. (13.171) can be ignored 
and only the saddle point in Eq. (13.151) contributes, corresponding to a principle of least 
action 6S = 0. In this way we recover the results in Section II C. We note that the canonical 
phase space method is completely equivalent to the Freidlin-Wentzel scheme for the extremal 
orbits. In fact inserting Eq. (13.11) for one variable, dx/dt = — F + p in Eq. (12.211) we obtain 
S* = (1/2) J dtp"^ which is the Freidlin-Wentzel action. The advantage of the phase space 
method is the introduction of the canonically conjugate momentum p, representing the noise, 
as an additional variable. This allows for a phase space representation of the numerical 
results obtained by a numerical optimization of the Freidlin-Wentzel action. In Fig. [6] we 
have in a xt plot depicted the paths in configuration space from an initial configuration Xi 
at time t = to a final configuration x at time t = T. We have also shown the extremal 
path which dominates the path integral in the limit A 0. 

IV. MINIMUM ACTION METHOD FOR THE KPZ EQUATION 

In this section we apply the minimum action method to the KPZ equation and set up 
the numerical scheme. For the KPZ equation the FW action has the form 



(4.1) 
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In order to find the optimal switching path from an initial configuration /linit(r) at time 
t = to a final configuration /ifin(r) at time T we minimize the action (14.11) subject to the 
constraints: 

/i(r,0) = /ii„it(r), Mr,T) = /ifi,(r). (4.2) 

We first discretize the action functional using finite differences, then minimize the discretized 
action functional using the limited memory Broyden-Fletcher-Goldfarb-Shanno (BFGS) 
method. The BFGS is an efficient quasi-Newton method for large scale optiminization 



problems 7^. It is an iterative method; at each iteration, it only requires the input of 
the action S and the associated gradient 6S/6h{r,t). The minimization is constrained by 
appropriate Dirichlet boundary conditions in space, h{r,t) = hig{r) for r on boundary B, 
and initial and final boundary conditions in time, h{r,0) = /iinit(i') and h{r,T) = /ifin(i") 

In the following we consider the ID case. We confine the system to a ID interval of size 
L and the switching path to the time interval T. In the space-time domain [0, L] x [0, T] we 
introduce a mesh with sizes Ax = L/I and At = T/J and define the grid points {xi,tj) 

Xi = iAx ii = 0, ■■■/, (4.3) 
t, = J At J = 0, ■ • ■ J. (4.4) 

The numerical approximation to h{xi, tj) is denoted by Hij. In order to simplify the expres- 
sion we introduce the momentum or the noise field 

, , dh d^h \fdhV ^ 
p(a;,t) = --^.— - +F, (4.5) 



dt dx"^ 2 \dx J 

and express the action in the form 



S{h) = 11 dt [ dxp^{x,t). (4.6) 
2 Jo Jo 



Using the trapezoidal rule to discretize the integral in space and the midpoint rule to compute 
the temporal integral we obtain 

/-I J 

AxAt 

2 



S{H) = lAxAtJ2J2Pi^ (4.7) 



i=i j=i 
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where the discretized version of the noise field is 




Hi+i,j + Hi_ij — 2Hij 





(4.8) 



For the discretized boundary condition we have 



Hoj^H,, Hij^H2 for j = 0, • • • J 

Hio = hinitixi), Hij = hfinixi) for i = 0, ■ ■ ■ /, 



(4.9) 



(4.10) 



where Hi and H2 denotes the boundary values. For an offset in the height profile we have 
Hi ^ H2. The BFGS method also requires the gradient of the action, whose discrete version 
is given by 



The numerical optimization is set up by choosing an initial pathway interpolating between 
the initial and final configurations /iinit and /ifin subject to the chosen boundary conditions. 
Provided that the initial pathway lies in the domain of the appropriate minimum of S the 
BFGS method then through successive steps finds the minimum action and outputs the 
weak noise pathway from hir^aix) to /ifin(a;) in a given transition time T. 




2(A,t)2 

-Pi+ij+i + Pi-i,j+i — 2-Pi,j+i 



2{Axy 




(4.11) 
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V. NUMERICAL RESULTS IN ID 



In this section we discuss various switching scenarios for the KPZ equation in ID. As 
parameter values we choose for the viscosity u = 1 and for the non-hnear growth parameter 
A = 2. These values yield the inverse length scale ko = 1. The parameter k is then given by 
k = a/F where F is the imposed drift. We, moreover, consider a system of size L = 1. 

We consider the switching scenario in ID from an initial state h{x, 0) = —Hq to a final 
state h{x,T) = Hq. This transition corresponds to adding a layer of thickness 2ho to the 
interface. At the boundaries x = and x = L we set h = 0, i.e.. Hi = H2 = 0. In order to 
match the initial profile /iinit(2;) = —fiQ to the boundary condition we use the cusp solutions 
in Eq. (I2.37p . h±{x) = ±(l//i:o) In | coshfca;| and set h{x,0) = h^^x) + hR{x), where 



1 

hL{x) = -—In 

Kq 

hR{x)= -^In 



cosh k{x — Xi] 



cosh k{x — Xi — 5) 
cosh k{x — X2) 



(5.1) 
(5.2) 



cosh k{x — X2 — 5) 

Setting xi ~ and X2 ^ L — 5 and choosing 5 = hokQ/2k the initial profile satisfies the 
boundary conditions and approach the interface value —ho in the bulk; note that the slope 
of the steps is given by 1/k. In our simulation we choose xi = 0.1, X2 = 0.8, 6 = 0.1. 
Likewise, the final configuration at time t = T is given by h{x,T) = — /il(x) — hii{x). To 
ensure a steep step corresponding to a short healing length we choose the drift F = 625 
corresponding to A; = 25. With this choice ho = 2k6/ko = 5. For the initial path, we use the 
linear interpolation between h{x,0) and h{x,T): h{x,t) = (1 — t/T)h{x,0) + {t/T)h{x,T). 
Finally, we have chosen a 200 x 200 set of xt grid points. 

In Figs. [3-[in]we show switching scenarios for the transition times T = 0.1, T = 0.03, 
T = 0.01, and T = 0.001. We depict both the height profiles h{x,t), the slope profiles 
u{x,t), and the associated noise profiles p{x,t). In Figs. [TTIfT^ we depict the associated 
squared noise field or action density in a space-time plot. 

The height profiles presented for the initial and final configurations /imit and /ign and for 
some characteristic intermediate times show that the transition or switch in time T is effec- 
tuated by the ballistic propagation of steps or facets across the system. The corresponding 
slope profiles demonstrate that the steps can be interpreted in terms of a gas of domain walls 
with opposite parity, i.e., right hand and left hand domain walls. The motion of a single 
step in h is thus associated with a pair of co-moving domain walls in u moving across the 
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Transition time T 


Switching action S 


0.100 


2.57 X 10^ 


0.030 


2.56 X 10^ 


0.010 


3.12 X 10^ 


0.001 


1.95 X 10^ 



TABLE I: The switching actions >S'(r) associated with the transition times T = 
0.100,0.030,0.010,0.001. 

system. The dependent noise field p is associated with the nucleation of domain walls. Since 
the right-hand domain wall is a solution of the deterministic Burgers equation it carries 
no dynamical attributes and the associated noise field vanishes, unlike the "noise-induced" 
left-hand domain wall which is associated with a noise field and carries an action. In Table I 
we show the actions associated with the transitions and in Figs. [15] the action as a function 
of the transition time for the various scenarios. 

VI. DISCUSSION AND INTERPRETATION 

In this section we interpret the numerical results in ID using the weak noise canonical 
phase method. As discussed in Sec. |TTl]the phase space method is completely equivalent to 
the minimum action method. 

A. Waiting time transition for T = 0.1 

In terms of the switching dynamics T = 0.1 corresponds to a long-time transition. In 
Fig. [7]we show snapshots of /i, m, and p at times t = 0.0, 0.05, 0.0875, 0.0925, 0.1; in Fig. [TT] 
we depict the squared noise field or space-time action density. In the initial stage of the 
transition, from t = to about t = 0.075, the constant height field makes a transition to a 
trough (convex cusp) compatible with the boundary conditions h = 0. This configuration 
corresponds to a static right hand domain wall in the slope u. After a long waiting time 
in this configuration (until about t = 0.075) domain walls in u nucleate at the boundaries 
and a pair of domain walls then move across the system from left to right. In the height 
field this mode corresponds to the motion of a facet or step. The trough in h is filled in and 
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eventually at time T the final configuration /ign is reached. The noise field associated with 
the waiting time configuration vanishes since it corresponds to a right hand domain wall. 
However, in the final stage of the transition the noise field develops corresponding to the 
nucleation of the left hand domain wall. 

This switching scenario is in accordance with the phase space interpretation generically 
represented in Fig. [2l For a long time transition the orbit comes close to the saddle point 
corresponding to p = 0. In the slope field this implies a configuration given by the right 
hand domain wall u = {k / ko) tanh kx yielding the cusp in Fig. 7a. After a long waiting 
time in the vicinity of the saddle point the orbit eventually wanders off along the stationary 
manifold towards the final configuration. This part of the orbit in phase phase associated 
with a finite noise field corresponds to the propagation of the step in h, associated with the 
domain wall pair in u. 

The action can also be estimated qualitatively. For a single left-hand domain wall the 
action is given by Eq. ( ]2.34p . Sdw = {8/3)h''^{k^/kl)T. Inserting u = 1, ko = 1, and k = 25 
we obtain Sdw = 41667 x T. However, owing to the waiting time only the last p ^ part 
of the orbit contributes to Sdw Estimating the effective transition time to be T ~ 0.05 we 
obtain an action of order Sdw ~ 2000 which should be compared with the numerical value 
from Table I, S'num = 2567. The discrepancy can be accounted for by the finite nucleation 
action at the boundaries and also the finite system size effect. 

B. Intermediate time transitions, T = 0.03 and T = 0.01 

In Figs Eland [9] we have depicted switching scenarios at transition times T = 0.03 and T = 
0.01 for the height, slope and noise. In Fig. [8]we show snapshots along the pathway at times 
t = 0.0,0.015,0.0188,0.0225,0.03 and in Fig. at times t = 0.0,0.0025,0.005,0.0075,0.01. 
In Figs. [12] and [13] we depict the squared noise field or space-time action density. Since 
the imposed transition time is shorter compared to the previous case the waiting time is 
shortened. The transition again is driven by the nucleation and subsequent propagation of 
domain walls. In the case T = 0.03 domain walls in u are nucleated at the edges and the 
pair propagates across the system with a positive velocity similar to the waiting time case. 
In the case T = 0.01 the shorter transition time favors the nucleation of a domain wall in 
u at the center. This domain wall subsequently breaks up into two pairs of domain wall 
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moving toward the edges. In the height profile this scenario corresponds to the nucleation 
of a tip which subsequently broadens to a plateau effectuating the transition. 

This switching scenario is again heuristically in agreement with the phase space inter- 
pretation in Fig. [21 For an intermediate time transition the orbit in phase space bends 
off towards the stationary finite p manifold at an earlier stage in order to effectuate the 
transition in the shorter time interval available. 

The action based on Eq. fl2.34p is again of the same order of magnitude as the numerical 
results listed in Table I. We note that the shorter transition time requires a larger domain 
wall velocity v oc ki, where ki is the charge of the particular domain wall. Since the action 
scales with kf this effect compensates in the action for the smaller T. For an infinite system 
the imposed drift F oc fc^ in the KPZ equation is related to the domain wall charges ki by the 
relationship k = ki. Due to finite size effects this relation cannot be used directly in the 
present context. However, we still conclude that the imposed k does not fix the individual 
charges. The domain wall amplitudes and thus velocities are determined by the transition 
scenario. 

C. Short time transition, T = 0.001 

In Fig. [in] we show the switching scenario for the transition times T = 0.001 for the 
height, slope and noise. In Fig. we show the snapshots along the pathway at times 
t = 0.0,0.00025,0.00075,0.001. In Fig. [14] we depict the squared noise field or space-time 
action density. In the short time regime it is more advantageous to nucleate multiple domain 
wall pairs in the slope field, corresponding to multiple steps of facets in the height field. 

D. Switching action 

In Fig. [15] we depict the action S{T) as a function of the transition time T for five 
transition scenarios. The circles correspond to the transition pathways we discussed earlier 
and shown in Fig. [TllTOlfor T = 0.01, 0.03, 0.01, 0.001; the remaining pathways (not shown) 
involve one nucleation at the center and one nucleation from the boundary. The plot clearly 
indicates that more domain wall pairs, yielding a lower action, are nucleated at shorter 
transition times. 
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This relationship can be accounted for by the following considerations. For a single 
domain wall pair propagating across the system the associated action is given by Si = 
Snuci + /T^T . Here S^^ci is the nucleation action associated with the left handed domain 
wall. The second term follows from Eq. fl2.34l) . where we note that the velocity v = L/T 
scales with the amplitude k; A is a. constant which we do not have to specify further. In the 
case of a transition effectuated by the nucleation and transition of two domain wall pairs 
we have correspondingly for the action, 5*2 = 2S'nuci + 2A{L/2T)^T, where we note that the 
domain wall pair only propagate half the distance. In the general case of n domain wall 
pairs we obtain the expression 

Sn = nSn^d + A——. (6.1) 

In Fig. [in] we have depicted >S'(T) versus T for different values of n, which shows that the 
multi domain wall transitions have lower action at shorter time. This result follows from the 
competition between the nucleation action and the action associated with the propagation 
and is in qualitative agreement with the numerical results shown in Fig. [151 



VII. TRANSITION PATHWAYS IN 2D 



In 2D the weak noise approach yields elementary spherically symmetric growth modes. In 
terms of the diffusive field w the diffusion equation fl2.3ip . V^w = k'^w, has the asymptotic 
growing solution w+ oc exp(A;r) for r ^ \jk giving rise to the height field /i+ = [kjki^T 
and the slope field u+ = {kjk^^Yjr. Likewise, the non-linear Schrodinger equation (12.321) . 
V^w = k'^w — k^w^, yields the decaying solution w_ oc exp(— A;r) and, correspondingly, 
/i_ = —{k/ko)r and u_ = —{k/kQ)v/r. The height modes correspond to a tip (upward cone) 
and a dip (downward cone) in the interface profile, whereas the slope modes are outward 
and inward pointing vector fields of constant magnitude k/k^, i.e., monopole fields. Like 
in the ID case the static growth modes can be boosted to a finite propagation velocity 
and one can construct a dynamic growth morphology in terms of a dilute gas of monopoles 
in the slope field with superimposed diffusive modes. In a charge language the positive 
monopoles are solutions of the noiseless Burgers equation and carry no action, whereas 
the negative monopoles carry an action S oc (z/^T/Zcq)/;;^. In order to model a pathway 
from /linit at t = to /ifin at t = T one assigns a gas of monopoles representing /linit- 
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With the appropriate assignment of the noise field corresponding to nucleation events this 
configuration will evolve in time to /ign- The total action associated with negative growth 
modes, using the WKB ansatz P oc exp{—S/A), then yields the transition probability for 
the kinetic pathway. Details of this procedure has been discussed at length in Ref. jssi and 
will not be reproduced here. 

The minimum action method is easily extended to higher dimension generalizing the 
procedure in Sec. IIVI Choosing the parameters z/ = 1, A = 2, and a 100 x 100 x 100 set 
of xyt grid points and matching the height profile to the boundary values h{r) = by a 
2D generalization of Eqs. (15. 1115. 21) . we have in Figs fT71 and fTSl depicted the 2D switching 
scenarios for the height field at transition times T = 0.02 and T = 0.002 from an initial 
plateau ai h = —5 to a final plateau at h = 5. In the case T = 0.02 a single peak in h 
is nucleated at the center of the plateau h = —5. The peak amplitude evolves in time and 
eventually flattens to the plateau at h = 5. In the case T = 0.002 the transition takes place 
subject to the nucleation of a regular pattern of growing cones in h which eventually broadens 
and merge together. Like in the ID case we note again that more peaks are nucleated at 
shorter transition times. 



VIII. SUMMARY AND CONCLUSION 



In the present paper we have applied the minimum action method based on the Freidlin- 
Wentzel scheme for rare events driven by weak noise to the KPZ equation for a growing 
interface. The KPZ equation is a non-gradient system and a characterization of kinetic 
pathways in a free energy landscape is not available, on the other hand, the pathways can 
be characterized as taking place in an action landscape. Correspondingly, the transition 
probabilities are characterized by the associated action of a specific pathway, unlike the free 
energy case for gradient systems where free energy considerations apply in the evaluation of 
the Arrhenius factor for the transition. 

The minimum action method basically identifies the kinetic pathway in the action land- 
scape by seeking a minimum of the action using an optimization technique. Once the mini- 
mum has been reached the method provides the kinetic pathway subject to given initial and 
final configurations combined with appropriate boundary conditions. We have conducted 
a detailed analysis of the ID case and find that the pathways can be characterized by the 
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nucleation and subsequent ballistic propagation of growth modes. These growth modes cor- 
respond to moving facets or steps in the KPZ height field and to moving domain walls in 
the slope field. We also find that the numerical results are in good qualitative agreement 
with the canonical phase analysis previously developed for the KPZ equation. We have, 
moreover, applied the minimum action method to the 2D case. 

In conclusion, we believe that the minimum action method provides a new tool in an- 
alyzing the kinetics of spatially extended or field theoretical non-gradient systems like the 
KPZ equation studied here. The method supplements previous scahng analysis of the KPZ 
equation in focussing on the pattern formation or many body aspects of kinetic transitions 
in the weak noise limit. 
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FIG. 1: DRG phase diagram for the KPZ equation to leading loop order in d — 2. We plot the 
effective coupling strength g = A^A/i/^ as a function of the dimension d. In d = 1 the DRG flow 
is towards a strong coupling fixed point with scaling exponents C = 1/2, z = 3/2. Above the lower 
critical dimension d = 2 there is a kinetic transition line, delimiting a rough phase from a smooth 
phase. On the phase line z = 2 and C = 0- The weak coupling smooth phase is characterized by 
z = 2 and = (2 — d)/2. Above d = 1 the scaling exponents in the strong coupling rough phase 
are unknown 
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FIG. 2: Phase space representation of the weak noise method. We show the zero energy submanifold 
p = corresponding to the transient behavior and the submanifold p = 2F determining the 
stationary distribution. The manifolds intersect in a hyperbolic saddle point (SP). The infinite 
waiting time at SP corresponds to the long time Markoff behavior. We depict a finite time orbit 
from Xi to X in transition time T and an infinite time orbit passing through the saddle point. 
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FIG. 3: We depict the static domain walls in the slope field corresponding to the solutions of the 
diffusion and non linear Schrodinger equations for the diffusive field . In a) we show the right hand 
domain wall. This domain wall carries vanishing action and is identical to the viscosity smoothed 
shock waves of the noiseless Burgers equation. In b) we show the noise induced left hand domain 
wall carrying a finite action. 
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FIG. 4: We depict a three domain wall growth morphology. In b) we show the slope field composed 
of two right hand propagating domain walls and a single propagating left hand domain wall. In 
a) we show the corresponding moving facets in the height profile. Note that the charges of the 
domain walls add up to zero implying a fiat interface at the edges, however, allowing for an offset 
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FIG. 5: In b) we show a co-moving two-domain wall configuration in the slope u. This pair mode 
corresponds to a moving step or facet in the height field h depicted in a) . The mode carries a finite 
action associated with the left hand domain wall. 
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FIG. 6: We show paths (dashed) from an initial configuration xi at time t = to a final configu- 
ration X at time t = T contributing to the path integral. We also depict the extremal path (solid) 
dominating the path integral in the limit of weak noise. 
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FIG. 7: We depict the transition scenario for transition time T = 0.1. In a) we show the waiting 
time configuration and the propagating step in /i, in b) the quasi static right hand domain wall 
and the propagating domain wall pair in n, and in c) the corresponding noise field associated with 
the propagating left hand domain wall (arbitrary units). 
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FIG. 8: We depict the transition scenario for transition time T = 0.03. In a) we show the 
propagating step in h, in b) the propagating domain wall pair in u, and in c) the corresponding 
noise field associated with the propagating left hand domain wall (arbitrary units). 
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FIG. 9: We depict the transition scenario for transition time T = 0.01. In a) we show the emerging 
plateau in h, in b) the left hand domain wah associated with the appearance of the peak in h and 
the propagating domain wah pairs emerging from the center in u, and in c) the corresponding noise 
field associated with the nucleation and subsequent propagation from the center (arbitrary units). 



37 



8 

H 4 
X 

o 

m 
-4 



"T T" 



"T T" 



"T T" 



a) 



height profile, T = 0.001 



1 

50 
150 



200 




"T T" 



"T r 



T T" 




-200 
x104 



LU 4 
CO 



slope profile, T = 0.001 



1 1 

. c) 


1 1 1 1 1 

noise profile, T = 0.001 


1 1 


1 1 1 1 1 1 1 1 1 











0.2 



0.4 0.6 

X 



0.8 



FIG. 10: We depict the transition scenario for transition time T = 0.001. In a) we show the 
propagation of the multiple steps or facets in h, in b) the associated domain wah pairs in u, and 
in c) the corresponding noise field associated with the nucleation and subsequent propagation of 
domain wahs (arbitrary units). 
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FIG. 11: We plot the squared noise field p{x,t)^ or action density as a function of x and t in the 
case T = 0.1. The plot show the waiting time aspects of the transition scenario (arbitrary units). 
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Action profile (T=0.03) 
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Action profile (T=0.01) 




FIG. 13: We plot the squared noise field p{x, t)^ or action density as a function of x and t in the 
case T = 0.01 (arbitrary units). 
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Action profile (T=0.001) 




FIG. 14: We plot the squared noise field p{x,t)'^ or action density as a function of x and t in the 
case T = 0.001 (arbitrary units). 
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FIG. 15: We depict the action S{T) as a function of the transition time T for five transition 
scenarios. The circles correspond to the transition pathways for T = 0.01, 0.03, 0.01, 0.001; the 
remaining pathways involve one nucleation at the center and one nucleation from the boundary. 
The plot shows that more domain wall pairs, yielding a lower action, are nucleated at shorter 
transition times. 
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FIG. 16: The action S{T) given by Eq. is plotted as a function of T for transition pathways 
involving up to 4 domain walls pairs in the slope field. The labeling indicates the number of 
domain wall pairs. The lowest action and thus the most probable transition is associated with an 
increasing number of domain wall pairs at shorter transition times (arbitrary units). 
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2D Height profile (T = 0.02) 







FIG. 17: We depict a 2D long time transition scenario for the height profile from an initial plateau 
at h = —5 to a final plateau at h = +5. The transition time is T = 0.02. The transition takes 
place subject to the nucleation of a single peak in h at the center. The peak eventually broadens 
as we approach the final configuration (arbitrary units). 
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2D Height profile (T = 0.002) 




FIG. 18: We depict a 2D short time transition scenario for the height profile from an initial plateau 
at h = —5 to a final plateau at h = +5. The transition time is T = 0.002. The transition takes 
place subject to a regular pattern of 9 nucleation zones. The peaks eventually broaden as we 
approach the final configuration (arbitrary units). 
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